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ABSTRACT 

The effect of rarefaction acceleration on the propagation dynamics and structure 
of relativistically hot jets is studied through relativistic hydrodynamic simulations. We 
emphasize the nonlinear interaction of rarefaction waves excited at the interface be- 
tween a cylindrical jet and the surrounding medium. From simplified one-dimensional 
models with radial jet structure, we find that a decrease in the relativistic pressure due 
to the interacting rarefaction waves in the central zone of the jet transiently yields a 
more powerful boost of the bulk jet than that expected from single rarefaction accelera- 
tion. This leads to a cyclic in-situ energy conversion between thermal and bulk kinetic 
energies which induces radial oscillating motion of the jet. The oscillation timescale is 
characterized by the initial pressure ratio of the jet to the ambient medium, and follows 
a simple scaling relation Tosdiiation oc (/jet,o/^amb,o)'^^- It is confirmed from extended 
two-dimensional simulations that this radial oscillating motion in the one-dimensional 
system manifests as modulation of the structure of the jet in a more realistic situation 
where a relativistically hot jet propagates through an ambient medium. It is found that 
when the ambient medium has a power law pressure distribution, the size of the recon- 
finement region along the propagation direction of the jet in the modulation structure 
A evolves according to a self-similar relation A oc where a is the power-law index 
of the pressure distribution. 

Subject headings: galaxies: jets — hydrodynamics — methods: numerical — rela- 
tivistic processes — shock waves 



1. Introduction 

Relativistic jets are coUimated bipolar outflows that have a velocity almost equal to light 
speed. They are ubiquitous among astrophysical systems consisting of a compact object sur- 
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rounded by an accretion dis k, e.g., active galactic nuclei (AGNs; Ferraril Il998|) . microquasars 



Meszaros 



(IMirabel & RodriguezI Il999|) and the central engine of gamma-ray bursts (GRBs; IPiran 



2004 



20061) . Although there are many works that try to determine the mechanism by which 



relativistic jets are accelerated and coUimated, the process is still not understood. 



Aloy & RezzoUal (120061) recently reported that relativistic jets can be powerfully boosted 
along the interface between the jet and ambient mediu m, if the jet has su fficiently large veloc- 
ity and specific enthalpy, and is overpressured (see also lAloy et al.l l2003l . 2005). A rarefaction 
wave excited at the interface converts the relativistic thermal energy of the plasma into kinetic en- 
ergy and yields an amplification of the Lorentz factor in the jet-ambient medium interface. This 
type of boost, which we label rarefaction acceleration, is not possible in Newtonian dynamics, but 
is an inherent process in relativistic hydrodynamics. 

The relat i vistic mag netohydrodynam i c effec t s on the rarefaction a cceleration were studied by 
Mizuno et all (l2008b and I Aloy & Mimical (120081) . iMizuno et all ilOOm showed that, especially in 
the case of a magnetic field perpendicular to the jet direction, the boost bec omes more powerful 
than that expected from pure hydrodynamic models. lAloy & Mimical (|2008b studied the radiative 
output of magnetized jets bounded by anomalous shear layers in which jets are boosted due to the 
rarefaction wave. It was indicated theoretically that the boosted Lorentz factor is well described, as 
a function of t he initial pressure ratio between the jet and the ambient m edium, by a simple scaling 
law derived by IZenitani et al.l (|2010l) (see also iKomissarov et al.ll2010|) . 



Rarefaction acceleration of jets is a process that is expected to occur under realistic condi- 
tions. Since a large amount of energy is stored in the jet at its launch site, the jet would be initially 
overpressured. Actually, the jet from AGNs often appear to be overpressured with relat ivistic ve- 
locity even as it propagates through the interstellar medium dBicknell & Begelmanll 19961) . In some 
cases it has two components: a faster, lighter compo nent surrounded by a slower, denser compo- 
nent dGiroletti et al.ll2004l: iMeliani & Keppensll2009[) . The jet that emerges from the central engine 
of GRBs is also expected to be hot and overpressured when it propagates inside the progenitor star. 
Furthermore, the fluid of the boundary layer of the jet would be boosted by rarefaction acceleration 
when the relati vistic jet breaks out on the surface of the progenitor star and the external pressure 
support drops (|Tchekhovskoy et al.ll2010h . The Lorentz factor of the jet can be boosted when it 
enters the low pressure region out side of the thick a ccretion torus which is a possible candidate of 
the central engine of short GRBs (|Aloy et al.ll2005() . 



In this paper, we study in detail how the rarefaction acceleration affects the propagation dy- 
namics of relativistically hot jets through one-dimensional (ID) and two-dimensional (2D) rela- 
tivistic hydrodynamic simulations. The nonlinear interaction of rarefaction waves excited at the 
interface between the jet and ambient medium is especially focused upon because it might have a 
potential impact on the boosting process and even alter the dynamics and structure of the jet. This 



-3- 



paper is organized as follows. In Section 2, the numerical results of our ID relativistic hydrody- 
namic simulations are presented. We perform 2D simulations in Section 3. Finally, we discuss and 
summarize our findings in Section 4 and 5. 



2. One-dimensional Simulations 



2.1. Numerical Models and Setup 



In order to investigate the interaction of rarefaction waves excited at the jet-ambient medium 
interface, we initially set a jet surround by ambient gas in the calculation domain, as is schemat- 
ically shown in Figure \T\ We solve the special relativistic hydrodynamic (RHD) equations in a 
one-dimensional axisymmetric cylindrical coordinate system. The fluid velocity has two compo- 
nents: the normal velocity v,. and tangential velocity Vj to the interface which separates the jet and 
ambient medium. Derivatives of the physical variables in the z-direction are assumed to be zero. 
We adopt, as the equation of state for the relativistic gas, the ideal gas law with a ratio of specific 
heats r = 4/3. The basic equations are then. 
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where 7 = 1/a/1 -(Vx/c)^-(v./cy is the Lorentz factor and h 
enthalpy. The other symbols have their usual meanings. 



(1) 
(2) 

(3) 

(4) 

c^ + rP/(r-l)p is the specific 



For our initial conditions, we assume a relativistically hot jet with larger pressure and smaller 
rest-mass density than the ambient medium. The initial density and pressure in the jet are chosen as 
Pjet.o = 0.1 and /jet,o = 1 , respectively. Those of the ambient medium are pamb,o = 1 and Pamb.o = 0.1. 
In addition, the jet velocity to the z-direction is relativistic Vjet,o = 0.99c, with a Lorentz factor of 
7jet,o ~ 7. The ambient medium does not move and the normal velocity v, is set to be zero initially 
in the calculation domain. In this paper, we set a simple model for the jet-ambient medium system 
in order to investigate the basic physics of the interaction of rarefaction waves excited at the jet- 
ambient medium interface. 

The normalization units in length, velocity, time, and energy density are chosen as the initial 
jet width Wjstfl, light speed c, light crossing time over the initial jet width W^et,o/c, and rest mass 
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energy density in the ambient medium pamb.oc^. The computational domain spans < r < 2. A 
uniform grid with a grid size Ar = 10"^ is adopted for our calculations. We use a reflecting bound- 
ary condition on the axis r = 0. The outer boundary of the grid uses the outflow (zero gradient) 
boundary condition. 



We use a relativistic HLLC scheme for hydrodynamics to solve the RHD equations (IMignone & Bodo 



2005h. The primitive variables ar e calculated from the conservative variables following the method 
of iMignone & McKinneyl (120071 ). We use a MUSCL-type interpolation method to attain second- 
order accuracy in space while the temporal accuracy obtains second-order by using Runge-Kutta 
time integration. 



2.2. Results of One-dimensional Simulations 

2.2.1. The Jet Boosting Mechanism ofAloy & Rezzola 2006 

Since the jet is initially hotter and has a higher pressure than the denser, colder ambient 
medium, three types of hydrodynamic wave are excited at the jet-ambient medium interface. A 
shock wave propagates outward through the ambient medium, and a rarefaction wave starts to 
travel toward the center of the jet. Behind the shock wave, there exists a contact discontinuity, that 
is an entropy wave, corresponding to the edge of the jet. The pressure decrease in the jet due to the 
rarefaction wave accelerates the gas in the jet-ambient medium interfac e in the tangential direct ion 



as a natural outcome of relativistic hydrodynamics, as was reported by lAloy & RezzoUal (|2006|) 



Following the analysis of IZenitani et al.l (120101) . by combining Equations ([3]) and (01) we can 
obtain 

7 ph\ TT+Vr- 



dt dr J c c dt 

dP 

~ -— C.-v^^c). (5) 

This indicates that a time-decreasing pressure is responsible for the acceleration of the gas in rel- 
ativistic hydrodynamics, unlike the non-relativistic case. The conversion from the thermal energy 
to bulk kinetic energy of the jet is then constrained by the relativistic Bernoulli equation which 
provides enthalpy conservation, 

7/1 ~ const. (6) 



The relativistic jet is accelerated by this mechanism, found by 
numerical model. 



Alov & RezzoUa 



(120061) . in our 
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2.2.2. Temporal Evolution of the System 

Figured shows the temporal evolution of our jet-ambient medium system. The color contour 
represents the spatial distribution of (a) the density p, (b) the pressure P, (c) the radial velocity v,-, 
(d) the tangential velocity v^., and (e) the Lorentz factor 7. Snapshots of the spatial distribution of 
the hydrodynamics variables at the time t=5, 27, 43 and 57 are illustrated in Figure |3l |4l |5]and 
|6l respectively. The dashed line in these figures denote the profiles of hydrodynamic variables at 
the time t = 0, 17, 33, and 50 for the better demonstration of the propagation of the rarefaction and 
shock waves. In the early evolutionary stage (0 < f < 5: phase (i)), described in section 3.3.1, the 
single rarefaction acceleration of the gas is observed in the jet-ambient medium interface. Figure |3] 
gives the spatial distribution of the jet-ambient medium system at ? = 0, 1 and 5. The Lorentz factor 
of the fluid at the jet-ambient interface is boosted to 12 due to the rarefaction acceleration when 
t=l and 5. One can find that the contact discontinuity moves outward, that is, the jet expands since 
the contact discontinuity corresponds to the edge of the jet. Note that the head of the rarefaction 
wave intersects the jet axis while the Lorentz factor of the fluid inside the jet reaches its maximum 
value at the tail of the rarefaction wave when t = 5. 

Subsequently at ? ^ 5 the rarefaction waves converge on the central region of the jet and an 
incident shock wave is excited at the tail of the rarefaction wave in the cylindrical jet, bringing 
a substantial change in the dynamics (phase (ii) in Figure |2l). The shock wave is excited at the 
tail of the rarefaction wave as a natural result of the cylindrical shock tube problem in which the 
characteristic lines intersect. The shock wave is located at r = 1 when ? = 17 as shown in Figure 
4. The gas pressure in the interacting region of rarefaction waves, confined by the surface of the 
shock wave, is then further reduced and becomes lower than that of the ambient gas (see Figure S]). 
Since, according to Equation ©, the thermal energy is converted to the bulk kinetic energy of the 
gas, the Lorentz factor of the gas in the interacting region is further boosted. The peak Lorentz 
factor of the gas inside the rarefaction region reaches 7 ~ 60 at f = 44 (see Figure |5]), which is a 
factor ~ 5 higher than that due to the single rarefaction acceleration at the jet- ambient medium 
interface. The boosted Lorentz factor of the fluid at the jet-ambient medium interface due to the 
single rarefaction acceleration is roughly 12 as shown in Figure 3, 4, 5, and 6. 

The interaction of rarefaction waves generates a strong inward pressure gradient behind the 
jet-ambient medium interface which acts to decelerate the radial expansion of the jet, turning ex- 
pansion of the jet into contraction around time t = 21 (see Figure |2l a), |2lc)). The contraction of 
the contact discontinuity results in converging flows inside the jet in the phase (iii). 

When t = 44, the inward shock waves proceeding the converging flow collide with each other 
at the center of the jet and propagate outward. The gas bounded by the shock is compressed 
and heated. Since according to Equation (|5]) a time-increasing pressure decelerates the tangential 
velocity of the jet, the Lorentz factor of the jet reduces and the specific enthalpy increases (see 
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Equation Q). Thus, in phase (iv), the Lorentz factor of the shock heated gas drops. The spatial 
distribution of the jet-ambient medium system at the end of the phase (iv) {t = 57) is demonstrated 
in Figure |6l 

When the shock wave encounters the contact discontinuity around t = 58, the system has 
almost returned to its initial state. The jet still has a sufficiently larger tangential velocity and 
specific enthalpy than the ambient medium, but the gas pressure in the jet, however, has become 
smaller than that of the initial state. This is the result of the energy conversion from thermal energy 
to bulk kinetic energy of the jet. Since the system is restored to a state that is almost the same 
as the initial conditions, the three types of hydrodynamic waves, an outward propagating shock 
wave, a contact discontinuity (the edge of the jet), and a converging rarefaction wave, appear at 
the jet-ambient medium interface. Therefore the contracting radial motion of the jet becomes an 
expanding motion. 



2. 2. 3. Oscillation of the Relativistic Jet 

After phase (iv), the system has returned to conditions similar to the initial conditions, and 
the system repeats the cycle of phases (i)-(iv) until the pressure of the jet becomes equal to that of 
the ambient gas. During the cycle, the radial motion of the jet oscillates between expansion and 
contraction, and the system undergoes successive exchanges between thermal energy and kinetic 
energy. Figure |7] shows the temporal evolution of (a) the jet width, (b) the maximum and average 
tangential velocity, (c) the maximum and average Lorentz factor, and (d) the average of the specific 
enthalpy in the jet. The oscillations of the Lorentz factor and the averaged specific enthalpy are in 
anti-phase in accordance with the relativistic Bernoulli relation Q. 

When the pressure inside the jet becomes almost equal to the ambient gas around time t = 
1000, the oscillation of the jet ends. The spatial distributions of the density, pressure, radial and 
tangential velocity and Lorentz factor at t = 2000 are depicted in Figure |8] At the quasi-steady 
state, there exists an accelerated region localized to the boundary layer of the jet, as is shown in 
Figure [8ld) and (e). The physical parameters in the interface of the jet are comparable to those 



expected from single rarefaction acceleration (|Aloy & Rezzollall2006|) . This is because the oscil 



lating motion during the relaxation stage, shown in Figure Ul has little influence on the physical 
parameters in the boundary layer. At the final quasi-steady state, the pressure in the jet is the almost 
same as that of the ambient medium. The initial relativistic thermal energy of the jet is converted 
into the bulk kinetic energy of the jet. 
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2.3. Scaling Law for the Oscillation Timescale 

As shown in the section 2.2.3, the initial non-equilibrium system evolves, through a transition 
stage, toward a quasi-steady state where a hydrostatic balance is established in the radial direction. 
During the transition stage, the oscillation timescale is almost constant while the oscillation am- 
plitude gradually decreases. In this section, we derive a scaling relation which can reproduce the 
oscillation timescale of our jet- ambient medium system. 

Figure Ha) indicates that, during the transition stage, the radial size of the jet oscillates around 
the jet width of the final quasi-steady state. We then approximate the typical oscillating width of 
the jet by the jet width of the saturated state illustrated in Figure [8] The typical oscillation time of 
the jet would be determined by the propagation time of the sound waves over the typical oscillating 
width of the jet. When using the physical parameters at the quasi-steady state, the oscillating time 
is evaluated as 

T = 7jet,sW^jet,s/Cs (7) 

in the laboratory frame. Here, the subscript s stands for the physical values in the final steady 
state. Since the thermal energy inside the jet is relativistic, the sound speed is well approximated 
by Cs ^ c/v^. 

In order to derive the typical oscillation time of the jet-ambient medium system, we need to 
estimate the average Lorentz factor in the jet 7jet,s in the quasi-steady state. Since the total energy 
in the jet is almost conserved during the transition phased, we can obtain the following relation 
neglecting the rest mass energy, which is smaller than the relativistic thermal energy in the jet, 

^jet,s^7/et,s'F'amb,0 = W^jet,0^7jet,0^^jet,0 • (8) 

Note that we replace the pressure in the jet in the steady state Pjet.s with that of the initial ambient 
medium Pamb.o because a pressure balance is established between inside and outside the jet in the 
final quasi-steady state. 

From Equations (|7]) and ([8]), we can give the scaling law for the oscillation time 




In Figure |9l we plot the oscillation time averaged over ten cycles for numerical runs with different 
initial pressure ratios. The solid line represents the analytic scaling we derived. This indicates that 
our numerical results are well captured by our simple scaling law shown in Equation 



'Our models indicate that the loss of the total energy is less than 0.4% in the transition phase 
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yiagnetic fields, which are expected to exist in relatiy i stic jets from AGNs and GRBs (iBlandford & Payne 



1982; IUchida& Shibatalll985l : IShibata & Uchidalll986l : lBlandfordll2000l : iLvutikov & Blandford 



20031), may alter the oscillation mechanism we have presented. In the poynting flux jet, the mag- 



netic pressure is dominant over the gas pressure. Since th e sum of the magn e tic pressure and 
the gas pressure contributes to the rarefaction acceleration (IMizuno et al.l l2008l : lAloy & Mimica 



20081 : iZenitani et al.ll2010l : iTchekhovskov et all 120 101 : iKomissarov et al.ll201ol) . it is expected that 



the magnetic pressure would play much the same role as the relativistic thermal pressure in the jet 
oscillation. 



3. Two-dimensional Simulations 

3.1. Numerical Models and Setup 

We investigate how the oscillating motion found in the ID model affects the propagation 
dynamics and the structure of the jet in a more realistic two dimensional axisymmetric system. 
Our axisymmetric simulation of the jet propagation has been carried out in cylindrical coordinates 
(r,z), where the z-axis coincides with the symmetric axis (see Figure [T). Relativistically hot flow 
is continuously injected into the ambient medium from the lower boundary of the computational 
domain, z = Z\ow The goveming equations we solved are 
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(11) 

(12) 
(13) 



where we use the same normalization units and symbols defined in section 2 in order to compare 
the results of ID and 2D calculations. The injection flow has the same hydrodynamic parameters 
as the ID model, that is Pjet,o = 0.1 , ^jeto = 1-0, Vr,o = 0, v,,o = 0.99c, and 7jet,o ~ 7. 

For the simulation presented here, we use simple power law models for the ambient medium. 
When assuming a polytropic atmosphere, the ambient medium has pressure and density distribu- 
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tions 



■Pamb - -Pamb.ol T7} ) (14) 

"jet.O / 

f \-3a/4 

Pamb Pamb.O 

(77^ , (15) 



jet.O 



where r = y/r^ + z^ is the spherical radius in the cylindrical coordinate system and Pamb.o = 0.1 and 
Pamb.o = 1- Note that the model with a = corresponds to a uniform ambient medium model. For 
simplicity we neglect the effect of the gravity and the ambient medium does not move initially in 
the calculation domain. We set the power law index of the ambient medium a as 0.0, 0.4, and 0.8 
in the following simulations. 

The calculation domain spans (lOOH^et.o x 1500H^et.o) in the (r x z)-plane which corresponds 
to a 1600 X 24000 grid. A uniform resolution of 8 numerical cells over the radius of the injection jet 
is used. At the lower boundary (ziow = 1) hydrodynamic variables are fixed inside the jet injection 
region (0 < r < 0.5) while the boundary conditions are reflective outside the jet injection region. 
An outflow boundary condition is imposed on the outer boundaries of the grid and the symmetry 
axis is reflecting. 



3.2. Analytic Estimation: The Size of Cusp-Shaped Boosted Region 



Our one-dimensional models suggest that the relativistically hot jet alternates between accel- 
eration and deceleration phases with the oscillation period derived in Equation <^ when it propa- 
gates through a uniform ambient medium. This is a consequence of the in-situ conversion between 
thermal and bulk kinetic energies inside the jet. When the relativistically hot jet is continuously 
injected into the ambient medium in more realistic multi-dimensional situations, we can expect 
that the radial oscillating motion that appeared in the one-dimensional system manifests as the pe- 
riodically modulated structure of the jet along the jet propagation direction. From the scaling law 
for the oscillation period of the jet obtained in section 2.3, we estimate here the typical size of the 
region where the relativistic jet is boosted. 

The space-time diagram of Figure |2] shows the boosted region of the jet has a cusp shape 
which corresponds to the interacting region of rarefaction waves. The multi-dimensional extension 
of this jet oscillation in the ID system wo uld provide the periodic forrnation of the cusp-shaped 



boosted region confined by oblique shocks (|Norman et al.lll982l : ISanderslll983l) along the propaga 



region are observed inside relativistic jets, that is called the reconfinemet shock (e.g.. 


Marti et al. 


1997; 


Gomez et al. 


1997 


Komissarov & Falle 


1997 


1998; 


Alov et al. 


2000a 


b; 


Agudo et al. 


2001; 



Zhang et alJ2003[l2004l:lMizuta et aljl2006l:lMizuta & AlovllOOglJperucho & M artillOOvUM orsonv et al 



2007 . 



2OI0I: Lazzati et alJbood : iMimica et alj|2009l : lNagakura et alJbol The cusp-shaped boosted 



region, which is result from the jet oscillation due to the interacting rarefaction waves, would ap- 
pear in the multi-dimensional system as the reconfinemet region found in these previous works. 

We can approximate the size of the cusp-shaped boosted region, that is the reconfinemet 
region, by the propagation distance of the jet fluid during the typical oscillation time r written in 
equation Since the fluid velocity of the relativistic jet is almost equal to the speed of light, the 
propagation distance A is given by 



A = cr = V37jet,oW^jet,0 



Het.O 



amb,0 



1/2 



(16) 



The size of the reconfinement region is proportional to the Lorentz factor of the injected jet and 
the square root of the initial pressure ratio between the jet and ambient medium. When the rel- 
ativistically hot jet is injected continuously into the uniform ambient medium, the cusp-shaped 
reconfinement region with the size A will be formed periodically inside the jet. The scaling law 
which relates the typical size of the reconfinement regi on to the pressure of the ambient medium 
is also analytically derived by Daly & Marscheii (|l988|) . We compare our scaling to it later in the 
discussion section. 

In a gravitationally bounded atmosphere, which is plausible in the central engine of relativistic 
jets, the ambient gas should be stratified and have a pressure distribution along the propagation 
direction of the jet. Since the radial oscillation of the relativistic jet is controlled by the pressure 
ratio between the jet and ambient medium, the size of the reconfinement region A evolves with 
time when the jet propagates through the stratified medium. Adopting the pressure profile for the 
stratified ambient medium described by equation ([14] ). the size of the reconfinement region along 
the z-axis is evaluated as 



A = ySTjet.oW^jet.O 



Pi 



jet,0 



amb,0 



1/2 



jet,0 



a/2 



(17) 



Since the fluid velocity inside the reconfinement region should be relativistic, we can provide the 
following relation by replacing z in equation (flTI) with ct. 



A oc e'^ 



(18) 



This implies tha t the cusp-shaped reconfinement region evolves self-similarly in the multi-dimension al 
system (see a lso Komissarov & Fall j 1998 ; Bromberg & Levinsonll2007 ; Nalewajko & Sikorall2009 ; 
Kohleretallboill) ^ 



- 11 - 



3.3. Self-Similar Evolution of Reconfinement Region 

3.3.1. Uniform Ambient Medium Model 

Figure [lOl shows the temporal evolution of the injected relativistically hot jet for a uniform 
ambient model (power-law index a = in equations (fT4l) and ([15])) in the early evolutionary stage. 
The color contour represents the spatial distribution of (a) the density p, (b) the pressure P, and (c) 
the Lorentz factor 7 when t = 100 (left), 200 (middle), and 300 (right) respectively. As expected 
from one-dimensional models, the cusp-shaped rarefaction regions confined by oblique shocks are 
formed inside the jet. In those regions, the fluid of the jet is accelerated due to the interaction 
of rarefaction waves and subsequent in-situ energy conversion from relativistic thermal energy to 
bulk kinetic energy (see Section 2.2). Since the rarefaction waves are repeatedly excited behind 
the jet head, the reconfinement regions are periodically formed along the jet-propagation direction 
and create the modulatied structure of the jet in the early evolutionary stage. 

Figure [TTI gives the spatial distribution of (a) the density p, (b) the pressure P and (c) the 
Lorentz factor 7 when t = 2000. As the jet propagates, the modulated structure of the jet has a loss 
in coherency except at the region near the injection point. This is because there are vortices that are 
episodically created by a K elvin-Helmholtz like instability that develops between the jet and the 



backflow near the jet head ([Mizuta & Aloyll2009h . Since the continuously induced vortices make 



the ambient medium inhomogeneous, the modulation structure of the jet becomes incoherent and 
then the size of the reconfinement region violently fluctuates, except the region near the injection 
point where the ambient inhomogeneity remains relatively weak. 

The spatial distribution of the Lorentz factor near the jet injection point is shown in Fig- 
ure [T2tb). In order to compare it to the result of the ID calculation, the temporal evolution of the 
Lorentz factor in the jet-ambient system of ID simulation is shown in Figure [III a). The size of 
the cusp-shaped rarefaction region in both cases seems to be almost the same, although the vertical 
axis is the time in Figure [T2l a). The typical size of the reconfinement region in the jet-propagation 
direction can be estimated theoretically from relation ([T6l) . With using the parameters 7jet,o — 7, 
Wjetfi - 1, and Pjetfl/Pm-nb.o = 10, wc Can obtain A ^ 40, indicating our scaling relation captures the 
result of 2D simulation. The width of the contact layer of 2D simulation is larger than that of the 
ID calculation due to the lower resolution of the simulation. 

The pressure on the z-axis in the ID and 2D simulations is shown in Figure [12 a) and (b), 
respectively. Note that the horizontal axis represents the time t in Figure \\3\ a) and (b) depicts the 
pressure in a ID cut along the z-axis of Figure [TTT b). In ID models, the oscillation of the jet is 
a transient phenomena and terminated when the pressure inside the jet falls to the same level as 
that of the ambient medium (See section 2.2.3). On the other hand, in 2D case, the pressure inside 
the jet continues to vary and a pressure balance between the jet and the ambient medium is not 
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achieved. The predicted size of the reconfinement region apparently seems to exist only near the 
injection point. 

Figure [14] shows the time-distance diagram for the pressure along the z-axis. The recon- 
finement region near the injection point does not move with time while the reconfinement region 
formed behind the jet head propagates along the jet direction. Though the modulation properties 
of the jet change due to the ambient medium inhomogeneity, the typical size of the reconfinement 
region is expected to be maintained at a constant value at any time as shown in Figure [141 

The temporal evolution of the Fourier spectrum of the pressure measured along the jet axis is 
given in Figure \T5\ Here the Fourier transformation P(k) is obtained from 

1 /"^Ji^' 



where Ljet is the length of the jet in the propagation direction. The solid, dashed and dotted curves 
indicate the spectra at the different phases t = 600, 1200 and 2000, respectively. Despite the dif- 
ferent evolutionary phase, all the spectra show a peak at around /:/27r ~ 1/40 which corresponds 
to the inverse of the typical size of the reconfinement region derived from equation ([T6] ) and we 
can find that there is a second harmonic in each spectrum. This suggests that the typical size of the 
reconfinement region inside the jet is essentially determined by the modulation caused by the in- 
teraction of rarefaction waves even though the ambient medium inhomogeneity causes fluctuation 
in jet structure. 



There is one main important difference between the power law models with a ^0 and the 
uniform ambient model (a = 0). For the models with decreasing pressure profile with distance 
from the launching point of the jet, the jet expands more than in the uniform ambient model. The 
jet is then further accelerated because the pressure inside the reconfinement region decreases more 
drastically and a larger amount of thermal energy is converted into bulk kinetic energy. 

Figure [16] shows the temporal evolution of the injected relativistically hot jet for the power 
law ambient medium model with a = 0.8. The meanings of the color contours and the columns are 
the same as those of Figure [TO] Note that the vertical and horizontal scale of Figure [l6]is 2.5 times 
larger than that of Figure [10] It is found that, though the inflow flux is the same for each model, the 
maximum Lorentz factor inside the reconfinement region reaches ^180 for the power-law model 
with case a = 0.8, while it is ^ 45 for the uniform ambient medium model (see Figure [TOlc) and 
Figure [Htb)). 
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3.3.2. Power-law Models 
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The decreasing ambient pressure results in a longer oscillation timescale as suggested in equa- 
tion The size of the reconfinement region in the power law models thus becomes larger than 
that in the uniform ambient medium model. Additionally, the results shown in Figure [T6] imply 
that the reconfinement region evolves self-similarly. 

In order to confirm the self-similarity, we investigate the time trajectory of the cusp of the 
reconfinement region on the jet axis, depicted by diamonds in Figure [l7] after the cusp-shaped 
region is formed. Overplotted on Figure [T7]is the curve denoting the scaling relation (fTSl ) which 
is obtained from the results of the ID models (see section 3.2). Panels (a) and (b) correspond to 
the cases with a = 0.4 and 0.8 respectively. These figures confirm that the reconfinement region 
evolves self-similarly once the cusp-shaped reconfinement region is formed. 



4. Discussion 



Aloy & RezzoUal (120061) discussed, for the first time, a new mechanism in relativistic hydro- 
dynamics that can act as a powerful booster in jets, that is rarefaction acceleration. This mechanism 
is purely hydrodynamical and operates when the jet-ambient system satisfies three physical con- 
ditions, that is, the jet should be hot and overpressured, and have a relativistic velocity to the 
jet-ambient medium interface. The hot jet stores an internal energy comparable to or larger than 
its rest mass energy. When the pressure of the jet is larger than that of the ambient medium, the 
rarefaction wave is excited at the jet- ambient medium interface and the relativistic thermal energy 
of the plasma can be converted into the bulk kinetic energy of the jet near the interface. An ex- 
pansion of the relativistically hot jet to the underpressured ambient medium is responsible for the 
rarefaction acceleration of the jet. 

The rarefaction acceleration should be commonly operated in the jet which satis fies these con- 
ditions, and which have been observed in many previous simulati ons as expected (e.g. jGomez et al 
19971 : 1 Aloy et al.ll2000bl : IScheck et al.l l2002l : Izhang et al.ll2003h . The reconfinement region con- 
fined by the oblique shock has been also formed in these previous simulations as a natural outcome 
of the rarefaction acceleration in the multidimensional system. However, the rarefaction accelera- 
tion can not be observed in the cold jet where the internal energy of t he fluid is less than the res t 
mass energy of the fluid even if it shows the reconfinement shock (e.g., Komissarov & Fallal 1997 ). 



Daly & Marscheii (|1988h explored the dynamics of relativistic, hot and overpressured jets on 



the basis of a simplified, quasilinear, hydrodynamic equations for adiabatic, steady, cylindrically 
symmetric, and irrotational flows. In that work, they have already found the variety of flow struc- 
tures with oscillating cross section or standing shocks. In addition, they have discussed scaling 
laws which relate the intrinsic properties of the jet to the pressure of the ambient medium. 
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In this sense, our finding is that there exists a close connection between the rarefaction ac- 
celeration mechanism discussed by Aloy & Rezzola and the flow structure with the reconfinement 
region and oscillating cross section quantitativel y by using ID and 2D hy drodynamic simulations. 
It is remarkable that simple ID model based on lAloy & RezzoUal (|2006|) 's mechanism can repro- 
duce complicated 2D structure of jets not only in rarefaction acceleration region but also in the 
reconfinement region so well. Since ID model is simple and easy to understand, these findings 
would be useful for more detailed study of jets. 

It should be finally emp hasized that the scaling la w we obtained in this work (equation (fT6l )) is 
different from that found by iDaly & Marscherl(|l988l) . This would be because the approximations 
used for deriving the scaling in their work is applicable only to the system with a small pressure 
difference between the jet and ambient medium. Our scaling law can capture the flow properties 
of the jet-ambient medium sys tem in the wider parameter range than that covered by the scaling 
law of balv & Marscheil (Il988l) . 



The simple scaling law for the self-similarity of the reconfinement region (IT8I I would pre- 
dict the evolution of the reconfinement region in the context of both AGN jets and GRBs. The 
reconfinement shocks due to the interacting rarefaction waves ma y be able to account for the con 



finernent of the relativisti c jet from AGN on a sub-parsec scale (|Junor et al .11 19991: iKovalev et al 



2003; lAcciari et al.ll2009|) . So me authors propose that the radio knots in the jet may correspond 



to reconfinement regions (e.g.. DDalv & Marscheiill988 



19961 : iKomissarov & Fallelll997l : IStawarz et al.ll2006l : 



Gomez et a 



Mimica et al 



19951 : iBickneU & Begelman 
2OO9I) . The scaling law ^ 



predicts that radio knots in the jet from AGNs evolve self-similarly with time depending on the 
power law index of the ambient pressure distribution when the relativistically hot and steady jet is 
formed near the central engine. 

Although the numerical work in this paper are based on simulations with a constant engine 
which powers the jet, the variability injected by the central engine also leads to inhomogeneity 
in relativistic jets ([Gomez et al.lll997l: iMimica et al.ll2009l: iMorsony et al.ll2010|) . This might be a 
promising mechanism for generating internal sh ocks in relativistic jets, accounting for the time 



variable emission prop erties of AGNs and GRBs ([Takahashi et al.ll2000l : |Piranll2004l: iMimica et al 



2005l : lMeszarosll2006l) . The relation for inhomogeneity in the jet between the variability injected by 
the central engine and the interaction of the relativistic jet with the ambient medium is also needed 
to be investigated in a more realistic situation. Actually, for the case of power-law ambient models, 
we need a larger calculation domain while retaining the same resolution to study the effect of the 
vortices developed around the jet head on the propagation dynamics of the jet. This is because 
the jet expands more rapidly and the modulation structure becomes larger as the power-law index 
a increases. However, these studies go beyond the scope of this paper, but will be studied in a 
separate paper. 
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5. Summary 

The nonlinear evolution of the interacting rarefaction waves excited at the cylindrical jet- 
ambient medium interface is studied through one-dimensional relativistic hydrodynamic simula- 
tions. It is found that an enhanced decrease in the relativistic pressure due to the interaction of 
rarefaction waves transiently yields a more powerful boost of the Lorentz factor of the bulk jet 
than that expected from a single rarefaction wave. The cyclic in-situ energy conversion between 
thermal energy and bulk kinetic energy is a natural relativistic outcome of the jet scenario studied 
and responsible for the radial oscillating motion of the jet. The oscillation timescale is character- 
ized by the initial pressure ratio of the jet to the ambient medium, and follows a simple scaling 

relation Tosdllation (/jet,o/-Pamb,o)'^^- 

It is confirmed from extended two-dimensional simulations that repeated excitation and con- 
vergence of rarefaction waves result in the alignment of the interacting regions of rarefaction 
waves, confined by oblique shocks, along the propagation direction of the jet when a relativis- 
tically hot jet propagates through an ambient medium. The evolution of the reconfinement region 
in which the fluid of the jet is powerfully boosted due to the interacting rarefaction waves has a 
self-similar property when the relativistic jet propagates through the ambient medium which has 
a power law pressure distribution. The evolution of the size of the reconfinement region along 
the jet direction A follows a simple scaling law A oc f"/^ where a is the power law index of the 
pressure distribution. Especially, in the uniform ambient medium model (a = 0) the typical size of 
the reconfinement region inside the jet is essentially determined by the modulation caused by the 
interaction of rarefaction waves. 
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A. Convergence Tests 



A.l. One-dimensional Riemann Problem with Transverse Velocity 

The ID relativistic Riemann problem with large tra nsverse velocity requires high resolution in 



order to numerically reso lve the complicated structure (|Zhang & MacFadyenll2006l : iMizuta et al 



20061 : iMizuno et al.ll2008h . We have tested our code on a ID (x-direction) Riemann problem with 
transverse velocity using uniform resolutions of Ax = 10"^ and 2 x 10"^ when the calculation do- 



main spans < X < I. We have used the initial conditions from lZhang & MacFadyenI (|2006|) with 
adiabatic index F = 5/3 as follows. 

Left state (0 < jc < 0.5): pL = 1, A. = 10^ v.vx = and v^x = 0.99c. 

Right state (0.5 < x < 1): pR = 1, Pr = 10"^ v,,r = and v. r = 0.99c. 

Figure [TSl denotes numerical and analytic soluti ons of this problem at t = 0.6. A nalytic solu- 
tions in the appendix are calculated with the code of iGiacomazzo & RezzoUal (|2006[) . Upper and 
lower panels correspond to the case with Ax= 10~^ and 2 x 10"^, respectively. In the low -resolution 
run, we can find the diffusion of the contact discontinuity. The position of the right-going shock 
front is not correctly captured. On the contrary, the left-going rarefaction wave is resolved with 
good accuracy. In the high-resolution run, the numerical solution captures the analytic solution 
well, although there is some undershoot in the transverse velocity Vj. at the contact discontinuity. 



A.2. One-dimensional Hydrodynamic Relativistic Jet Model 

We have also performed another convergence test using our jet model parameters in ID (x- 
direction) cartesian coordinate using a uniform high resolution of Ax = 2 x 10"^ and the resolution 
of Ax = 10~^ adopted in section 2. We set the initial conditions as follows. 

Left state (0 < jc < 0.5): pL = 0.1, Pl = 1, Vr,L = and v^,l = 0.99c. 

Right state (0.5 < x < 1): Pr = 1, Pr = 0. 1, Vv,r = and v^-.r = 0. 

Figure [19] shows the numerical and analytic solutions of our jet model at ? = 0.1 in the case 
with Ax = 2 X 10"^. The numerical solution resolves the hydrodynamic structure of the analytic 
solution very well. In the case with Ax = 10"^, the position of the right-propagating shock front is 
not correctly captured while the numerical solution is able to capture the hydrodynamic structure of 
the rarefaction wave and the contact discontinuity as shown in Figure |20l In this paper, we focus on 
the interaction of rarefaction waves excited at the jet-ambient medium interface. Since the shock 
wave propagates outward from the jet-ambient medium interface, the interaction of rarefaction 
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waves inside the jet is not affected by the outward-going shock wave. 

In our calculations in section 2, shock waves are excited inside the cylindrical jet. In order 
to investigate the effect of these shock waves, we have compared results of the calculation we 
performed in section 2 to those of the high-resolution calculation, in which a uniform grid with the 
grid size Ar = 2 x 10"^ is adopted in the cylindrical coordinate. Figure |2T]illustrates the temporal 
evolution of (a) the jet width, (b) the maximum and (c) the average Lorentz factor in the jet in both 
numerical runs. Note that horizontal axis corresponds to the time for the duration < t < 100 
demonstrated in Figure[7] The jet width in the case with Ar = 10"^ is up to a factor 1 .07 larger than 
in the high-resolution case because of numerical diffusion of the contact discontinuity. This leads 
to a 2 percent decrease in the average Lorentz factor inside the jet compared to the high-resolution 
case. We can find, however, only small differences of the maximum Lorentz factor in the jet 
between these numerical runs. Therefore, our choice of resolution does not impact significantly on 
our main results. 
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Fig. 1 . — The numerical setting for our ID simulations. A relativistically hot jet region is bounded 
by ambient gas. Derivatives of the physical variables in the z-direction and 6'-direction are assumed 
to be zero. We calculate the evolution of the jet only in the radial direction. 



-21 - 



80 



(a) 



Dcusitv 



40 - 



-2 



10 


80 
70 - 



1 1 


1 


(iv) 

contact 

tiiscoiitiiLuity 


contact 

diiscoiitimiity 














W 1 1 1 



log p 
n 



Pressure 



logF 



- -2 




10 12 

r 

Radial Velocity iv/c 



-2 -10 1 
r 

(d) Tangential Velocity 



10 




T 



T 



T 



60 |— rfiock '■^ 

I shock s 1 lodk 

40 — ""^^ 

30 - /' 
I 

20 - : 



sbuck 





Cl.l 


80 






70 






so 






50 




0.0 


■« 40 






30 






20 






10 


1 


-0.1 






-2 



12 -2 
(e) Lorentz Factor 





1.00 



0.99 



- 0.98 



0.97 



- 0.96 



L-l 0.95 



Fig. 2. — Temporal evolution of the jet- ambient medium system: (a) the density, (b) the pressure, 
(c) the radial velocity, (d) the tangential velocity and (e) the Lorentz factor. In panel (c), the 
rarefaction region is enclosed by dashed-lines. The tangential velocity of the ambient medium is 
zero in the panel (d). (A color version of this figure is available in the online journal.) 
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Fig. 3. — The spatial distribution of the jet-ambient medium system at ? = (dashed lines), t = I 
(solid lines) and t = 5 (dash-dotted lines): (a) the density, (b) the pressure, (c) the radial velocity, 
(d) the tangential velocity and (d) the Lorentz factor. 
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Same as Figure |3l but at ? = 17 (dashed lines) and 27 (solid lines). 
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Fig. 5. — Same as Figure |3l but at ? = 33 (dashed lines) and 43 (solid lines). 
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Fig. 6. — Same as Figure |3l but at ? = 50 (dashed lines) and 57 (solid lines). 
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Fig. 7. — The temporal evolution of (a) the jet width, (b) the maximum and average of the tangen- 
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(d) the average of the specific enthalpy in the jet. The initial Lorentz factor inside the jet is roughly 
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Fig. 8. — The spatial distribution of the initial non-equilibrium and final quasi-steady states of the 
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Fig. 10. — Time evolution of the injected relativistic jet into the uniform ambient medium: (a) the 
density, (b) the pressure and (c) the Lorentz factor. The left, middle, and right columns correspond 
to t = 100, 200, and 300, respectively. (A color version of this figure is available in the online 
journal.) 
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Snapshots of the spatial distribution of (a) the density, (b) the pressure and (c) the 
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Fig. 12. — Left panel represents the time evolution of the Lorentz factor in the jet-ambient medium 
system of ID simulation. Right panel demonstrates the spatial distribution of the Lorentz factor at 
t = 2000 in the uniform ambient model (a = 0) of 2D simulation. (A color version of this figure is 
available in the online journal.) 
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Fig. 13. — Panel (a): temporal evolution of the pressure at the z-axis in ID simulation. Panel 
(b): spatial distribution of the pressure along the z-axis at t = 2000 in the uniform ambient model 
(a = 0) of 2D simulation. 
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Fig. 14. — The time-distance diagram for the pressure on the z-axis in the uniform ambient model 
(a = 0) of 2D simulation. (A color version of this figure is available in the online joumal.) 




Fig. 15. — Temporal evolution of Fourier spectrum of the pressure measured along the jet axis 
P{k) =-j^ A-axis(z) exp"'''^ dz, where Ljet is the length of the jet in the propagation direction. The 
solid, dashed and dotted curves correspond to the case t = 600, 1200 and 2000, respectively. (A 
color version of this figure is available in the online joumal.) 
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Fig. 16. — Time evolution of the injected relativistic jet in the power-law pressure distribution 
model where the power-law index a = 0.8: (a) the density, (b) the pressure and (c) the Lorentz 
factor. The left, middle, and right columns correspond to t = 500, 1000, and 2000, respectively. (A 
color version of this figure is available in the online journal.) 



-35- 





10 I — I 

10 100 1000 10000 

t 

Fig. 17. — Time trajectory of the position of the cusp of the first reconfinement region on the 
z-axis is depicted by diamonds. Solid lines demonstrate the scaling law A oc we derived in 
section 3.2 where A and a are the length of the cusp-shaped reconfinement region along the z-axis 
and the index of the power-law pressure distribution of the ambient medium, respectively. Panel (a) 
and (b) correspond to the case with a = 0.4 and a = 0.8. 
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Fig. 18. — Solutions of ID Riemann problem with large transverse velocity at t = 0.6. Upper 
and lower panels correspond to the high-resolution (Ax = 2 x 10"^) and low-resolution (Ax = 10"^) 
cases, respectively. Diamonds and crosses represent the numerical solution while solid lines are 
analytic solution. (A color version of this figure is available in the online joumal.) 
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Fig. 19. — Solutions of Riemann problem of our jet model in ID (x-direction) cartesian coordinate 
with the resolution of = 2 x 10"^. Diamonds and solid lines illustrate the numerical and analytic 
solutions, respectively. (A color version of this figure is available in the online journal.) 
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Fig. 20. — Same as Figure [19] but the case with A;c = 10 ^. (A color version of this figure is 
available in the online journal.) 
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Fig. 21. — Temporal evolution of (a) the jet width, (b) the maximum and (c) the average Lorentz 
factor in the jet. Diamonds and crosses correspond to the resolution of Ax = 10"^ and Ax = 2 x 10"^, 
respectively. (A color version of this figure is available in the online joumal.) 



